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ABSTRACT 

The gravitational lens system CLASS B2108+213 has two lensed images separated by 4.56 
arcsec. Such a wide image separation suggests that the lens is either a massive galaxy, or is 
composed of a group of galaxies. To investigate the structure of the lensing potential we have 
carried out new high resolution imaging of the two lensed images at 1 .7 GHz with the Very 
Long Baseline Array (VLBA) and at 5 GHz with global Very Long Baseline Interferome- 
try (VLBI). Compact and extended emission is detected from the two lensed images, which 
provides additional constraints to the lensing mass model. We find that the data are consis- 
tent with either a single lensing galaxy, or a two galaxy lens model that takes account of a 
nearby companion to the main lensing galaxy within the Einstein radius of the system. How- 
ever, for an ensemble of global power-law mass models, those with density profiles steeper 
than isothermal are a better fit. The best-fitting profile for a single spherical mass model has 
a slope of 7 = 2.45^0 ^ . The system also has a third radio component which is coincident 
with the main lensing galaxy. This component is detected at milli-arcsecond scales for the 
first time by the 1.7 GHz VLBA and 5 GHz global VLBI imaging. However, the third radio 
component is found not to be consistent with a core lensed image because the radio spectrum 
differs from the two lensed images, and its flux-density is too high when compared to what is 
expected from simple mass models with a variable power-law density profile and/or a reason- 
able core radius. Furthermore, 1.4 GHz imaging of the system with the Multi-Element Radio 
Link Interferometric Network (MERLIN) finds extended lobe emission on either side of the 
main lensing galaxy. Therefore, the radio emission from the third radio component is almost 
certainly from an AGN within the main lensing galaxy, which is classified as an FRI type 
radio source. 

Key words: gravitational lensing - quasars: individual: CLASS B2108+213 - cosmology: 
observations 



1 INTRODUCTION 

Gravitational lensing is caused by the deflection of light from a 
distant object by matter along the line-of-sight to the observer. 
If the surface mass density of an intervening matter distribu- 
tion is above some critical density then multiple images of the 
background source are formed. The positions and flux-densities 
of these lensed images can provide a powerful probe of the en- 
closed mass and shape of the lensing matter distribution (e.g. 
ISchneider. Ehlers & Falcolll992l) . If this mass constraint is com- 
bined with an additional estimate of the mass at another radius, 
then the average matter density slope [7, where the matter density 
p(r) oc r -7 ] of the combine d luminous and dark matter distri- 
bution can be established (e.g. ISand et ai]|2005l : iKoopmans et al.l 
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2006; GavazzTetaD[2007). For galaxy-scale haloes in the mass 
range lO^-lO 12 Mq, the results from combined lensing and stel- 
lar kinematic data suggest that the inner 0.3-5 effective radii 
of galaxies are consistent with an isothermal mass distributi on 
(7 = 2.01+g g^ lTreu & Koopmansll2004lKoopmans et alfeOOot) . 

There are now over 100 examples of galaxy-scale lens- 
ing known. These systems have been fo und from systematic 
searches around po t ential lens galaxies (e.g . iFassnacht et alj[2004l : 
iBolton et al.l 120061 : iMoustakas et al.l |2007f) and fro m large sur- 



veys of the lensed sou r ce parent populatio n (e.g. Bahcall et al.l 
1 19921 ; lKingetai]|l999l : iBrowne et al) 120031 : iMvers et all |2003|) . 
Serendipitous discoveries of gravitational lensing from deep high- 
resolu tion optical and infra red imaging have also been made 
(e.g. IFassnacht et al] l2006bh . Typically, the image separations 
for these gravitational lens systems are between 0.5-1.5 arc- 
sec, with only a handful having image separations > 3 arc- 
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sec dLawrence et alj 1 1984 Iwisotzki et alj 1 19931; iMufioz et all 
200 ll; ISluse et alj 120031; lOguri et al.l |2004 iMcKean et al.1 l20ol: 



Bolton et alj|200q ; llnada et alj|2006l) . Since the enclosed mass of 
a lens system is proportional to the square of the image sepa- 
ration ( i.e. Me oc Ad 2 for a circularly symmetric mass distri- 
bution; iKochanekl [l99lh . such wide image separation lens sys- 
tems (i.e. > 3 arcsec) may be due to haloes which are an order 
of magnitude more massive than those of typical lens galaxies. 
As such, wide image separation lens systems could be used to 
probe the matter distribution at the top end of the mass function 
for galaxy-scale structures. Alternatively, the wide image separa- 
tion could be due to the lens galaxy being in an over-dense en- 
vironment, for example, in a group or cluster of galaxies. Recent 
imaging and spectroscopic surveys of the local environments of 
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; Auger et al. 2007;]). However, the total 



contribution of the group or cluster to the lensing mass distribution 
is thought to be no more than ~5 per cent for systems with image 
separa tions of ~ 1 arcsecond dMomcheva et al.|[2006l ; lAuger et al.1 
l2007ah . Studies of gravitational lens systems with larger image sep- 
arations may show an enhanced lensing mass distribution which 
has been boosted by the environment. An extreme example of this 
is cluster lensing w hich can produce image separations m uch larger 
than 10 arcsec (e.g. lQguri et al. 2004; llnada et alj2006h . 

The gravitational lens system CLASS B2108 +213 has two 
radio-loud lensed images separated by 4.56 arcsec jMcKean et al.1 
l2005h . This image separation immediately identifies B2108+213 
as an excellent opportunity to study a mass regime between 
the typical galaxy and cluster-scales. The system was discov- 
ered from 8.46 GHz imaging with the Very Large Array (VLA) 
during the Cosmic Lens All-Sky Survey (CLASS; iMvers et al] 
l2003l ; lBrowneetall2 003 ) . Follow-up observations with the Multi- 
Element Radio Link Interferometric Network (MERLIN) at 5 GHz 
found the radio spectra of the two lensed images (A and B) to be 
similar (a ~ 0.15 between 5 and 8.46 GHz, where S v oc p a ), 
with a flux ratio of Sb/Sa ~ 0.5. High resolution radio imaging 
with the Very Long Baseline Array (VLBA) at 5 GHz found the 
surface brightness of images A and B to be consistent with gravita- 
tional lensing. Not surprisingly, optical and infrared imaging found 
a massive elliptical galaxy (Gl) at the expected position of the lens. 
Furthermore, a companion galaxy (G2) was also found within the 
Einstein radiu s . A m ass model for the system was proposed by 
McK ean et ail (2005) which incorporated both Gl and G2 as sin- 
gular isothermal spheres with an external shear component. This 
model successfully reproduced the observed positions and flux ra- 
tios of the lensed images and required only a small contribution 
from the field (the external shear was ~2 per cent). However, given 
the limited number of observational constraints that were provided 
from the two lensed images, this model has no degrees of freedom. 
A spectroscopic survey of the environment around the main lensing 
galaxy has found at least 40 group/cluster galaxies at the same red- 
shift as Gl (z = 0.365). Also, the central stellar velocity dispersion 
of Gl was found to be 360 kms -1 (McKean et al. in preparation). 
These new data confirm that Gl is a massive galaxy in a dense 
group/cluster environment. 

Intriguingly, the radio imaging of B2108+213 with the VLA 
and MERLIN also detected emission from a third radio compo- 
nent (C) which is coincident with the lensing galaxy Gl. Lens 
theory predicts that an extended mass distribution should produce 
an ad ditional lensed image near the centre of the lens potential 
(e.g. iRusin et alj|200ll ; lKeetorJl2003h . The detection of such core 



lensed images is extremely rare for galaxy-scale systems because 
the image magnification tends to zero as the inner density profile 
approaches isothermal. However, the non-detection of core lensed 
images can place a very strong lower bound to the density pro- 
file of the lensing mass distribution (e.g.|Rusin et al.ll200ll ; iKeetorJ 
2003; Bovce et al . 2006; Zh ang et alj2007l) . There is currently only 
one accepted detection of a core lensed image, PMN J1632— 0033, 
whose lensing ga l axy has a power law slope of 7 ~ 1 .91 ±0.02 
JWinn et alj|2002l ; lwinn. Rusin & Kochanek| l2003. 2004). The na- 
ture of the third radio component of B2108+213 is not entirely 
clear. The radio spectrum of component C from 5 to 8.46 GHz 
is consistent with that of the lensed images, A and B. However, 
the flux-density of component C appears to be too large to be the 
core lensed image; the flux ratio is Sc/Sa ~ 0.1. Furthermore, 
component C was not detected with the VLBA which suggests that 
the surface brightnes s is not consistent with a core lensed image 
jMcKean et all2005l) . 

In this paper, we present new high-resolution radio imaging of 
B2108+213 with the twin aims of determining the nature of the 
third radio component and finding additional observational con- 
straints for the lensing mass model from the two lensed images. 
In Section [2] we show 1.4 GHz imaging of the system with MER- 
LIN which suggests that the emission from component C is from an 
AGN within the lensing galaxy Gl. High resolution imaging of the 
lensed images at 1.7 GHz with the VLBA and at 5 GHz with global 
VLBI (Very Long Baseline Interferometry) are also presented in 
Section [2] These new data show extended jet emission from both 
lensed images on mas-scales. We test gravitational lens mass mod- 
els for the B2108+213 system in Section|3] We discuss our results 
in Section|4]and present our conclusions in Section[5] 

Throughout, we assume a Qm — 0.3 and Qa = 0.7 flat- 
universe, with a Hubble constant of Ho = 70 km s _1 Mpc -1 . 



2 RADIO OBSERVATIONS 

In this section, radio imaging of B2108+213 with MERLIN, the 
VLBA and global VLBI are presented. 



2.1 MERLIN 1.4 GHz observations 

B2108+213 was observed with MERLIN at 1.4 GHz in two runs. 
The first set of observations were carried out for 14 h on 2005 
March 19 and 20 with all of the MERLIN antennas except for the 
Lovell telescope. The second observing run, which included all of 
the antennas, lasted for 14 and 8 h on 2005 April 15 and 16, re- 
spectively. 3C286 and OQ208 were used as the flux-density and 
polarization calibrators. A switching cycle of ~2 and ~6 min was 
used between the phase calibrator (B2103+213) and the lens sys- 
tem. The data were taken in a single IF and divided into 15 channels 
of 1 MHz width each. The data were taken with both right and left 
hand circular polarizations. The MERLIN pipeline was used for 
initial calibration and editing of the data within AIPsQ. 

The MERLIN imaging detected a combination of compact and 
very low surface brightness extended emission. The total intensity 
image was initially cleaned with the AIPS task IMAGR to subtract 
the higher surface brightness compact radio components (A, B and 
C). The residual image, which contained the low surface bright- 
ness emission, was then deconvolved with the maximum entropy 
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Figure 1. HST /-band image (F814W) of B2108+213 in grey scale overlaid 
with the contours at 1 .4 GHz using MERLIN. The lensed images A and B 
are coincident in the optical and the radio. The core component C coincident 
with galaxy Gl shows low surface brightness emission extending on either 
side in the radio. Galaxy G2 does not show any counterpart in the radio. 
The contours are (-3, 3, 6, 12, 24, 48, 96, 192)x0.03 mjy beam -1 (rms in 
the map). North is up and east is to the left. 



routine VTESS with a starting model of the central part of the resid- 
ual image smoothed with a circular Gaussian of full width at half 
maximum (FWHM) of 2 arcsec. Finally, the high surface bright- 
ness features were restored with the AIPS task RSTOR into the low 
surface brightness map and the combined total intensity map was 
produced. 

Fig. \H shows the total intensity map made from the combined 
datasets. Components A, B and C in this image are restored with 
a 0.252x0.165 arcsec 2 beam. For the first time we find extended 
emission on either side of component C spreading over an area 
of ~10x2 arcsec 2 . It appears as though the extended jet emis- 
sion originates from component C. The positions and the flux den- 
sities obtained from fitting Gaussian model components to A, B 
and C with the task JMFIT in AIPS are given in Table [T] The to- 
tal flux-density for the extended structure, measured by integrating 
over the region within the 3-cr limit, is ~ 70 mjy. Earlier mea- 
surements from the NVSS0 catalogue I Condo n et al.l 1 1998h sug- 
gest a flux density of ~ 53 mJy for the system. The map of lin- 
ear polarized emission shows only image A, with a polarized flux 
density of 0.6±0.1 mJy. The rms noise of the polarized image is 
~0.1 mjy beam -1 . Therefore, these observations were not sensi- 
tive enough to detect any polarized emission from image B, assum- 
ing that the polarized flux ratio between A and B is same as the flux 
ratio. 

2.2 VLBA 1.7 GHz observations 

B2108+213 was observed on 2002 June 19 with the VLBA at 
1.7 GHz. The aim of this observation was to image any low sur- 



Table 1. The flux densities and positions of A, B and C from fitting Gaus- 
sian model components to the MERLIN 1 .4 GHz map. 

Comp. RA Dec S peak S to tal 

(mas) (mas) (mJybeam -1 ) (mJy) 

A 0.0±0.2 0.0±0.2 16.3±0.8 16.4±0.8 

B 2141.3±0.4 -4026.7±0.4 7.1±0.5 7.7±0.4 

C 1434.4±6.0 -2915.6±6.0 0.5±0.1 0.6±0.1 



face brightness extended jet emission from the lensed images and 
to detect component C on mas-scales. The observing run lasted for 
7.5 h in total, with ~5.2 h spent observing B2108+213. The cal- 
ibrator B2103+213 was used for phase-referencing, with a 5 min 
time cycle between the lens system (3.5 min) and the calibrator 
(1.5 min). The data were taken in the left-hand circular polariza- 
tion through 4 IFs, each with 8 MHz bandwidth. The aggregate bit 
rate was 128 Mbs -1 with 2 bit sampling. All of the VLBA an- 
tennas were available with the exception of Los Alamos. The data 
were correlated at the VLBA correlator where each IF was divided 
into 16x0.5 MHz channels and averaged over 2-s time intervals. 
For wide-field imaging with interferometers, bandwidth smearing 
of the uv-data can cause radial smearing which reduces the peak 
flux density of a component far away from the phase centre. To 
limit the effect of bandwidth smearing, we correlated the data at 
two positions centred on images A and B. Note that for the set-up 
used here bandwidth smearing effects are significant at a radius 
of > 1 1 arcsec, which is well outside the maximum lensed image 
separation (4.56 arcsec). 

Both the A and B datasets were reduced individually using 
AIPS. The data were amplitude calibrated from the system tem- 
perature and gain curve values for each antenna, and corrected for 
the change in the parallactic angle. The data for B2103+213 were 
phase and amplitude self-calibrated and the solutions were applied 
to the lens system B2108+213. The mapping of B2108+213 was 
carried out using IMAGR, without further frequency or time averag- 
ing. We used an iterative process of cleaning and phase-only self- 
calibration, using a 9 min solution interval, to map the lens system. 

In the upper panel of Fig. [2] we show the naturally weighted 
maps of lensed images A and B, and of radio component C from 
dataset B. Image A shows a compact core and a radio-jet extending 
to the south-west which was fitted by a three component Gaussian 
model. Previous VLBA imaging at 5 GHz bv lMcKean et alj d2005h 
found the jet to extend over 10 mas in scale, whereas the new deeper 
1 .7 GHz map shows more emission extending beyond 20 mas from 
the radio core. Image B shows a single core component with a hint 
of jet emission to the north-east. Image B is also fainter than image 
A, which is consistent with gravitational lensing where the surface 
brightness of the lensed images is conserved. A single component 
Gaussian fit successfully modelled the emission from image B. The 
radio component C was detected for the first time at mas-scales. A 
single component Gaussian fit shows that the emission from C is 
unresolved. The results of fitting Gaussian model components us- 
ing the AIPS task JMFIT are listed in Table[2] Note that the extended 
emission detected in the MERLIN imaging has been resolved out 
here. 
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Figure 2. The top panel shows VLBA 1.7 GHz maps of B2108+213 restored with a beam size of 10.9x6.1 mas 2 and position angle of 12.18°. The bottom 
panel shows Global VLBI 5 GHz maps restored with a beam size of 3.19x0.97 mas 2 and position angle of —6.17°. Component A shows extended structure 
to the south-west direction at both frequencies. Component B shows north-east extension, which is better resolved at 5 GHz. Component C which is coincident 
with the main lensing galaxy, is compact at both frequencies. The contours for all of the maps are (—3, 3, 6, 12, 24, 48) X 

& map given in Table |2] for the 

1.7 GHz dataset and is 0.05 mJy beam -1 for the 5 GHz dataset. North is up and east is left. Grey-scales are in mJy beam -1 , except for the 1.7 GHz image of 
component C which is in /^Jy beam - 1 . 



Table 2. Positions, flux densities and rms noise of the components at 1.7 GHz. 



Comp. 


RA 


Dec 


Speak 


Stotal 


G map 




(mas) 


(mas) 


(mJy beam -1 ) 


(mJy) 


(mJy beam - 1 ) 


Al 


O.OrbO.l 


0.0±0.1 


3.0±0.2 


4.6±0.2 


0.070 


A2 


3.0±0.2 


-0.9±0.2 


2.3±0.1 


2.4±0.1 


0.070 


A3 


-12.2±0.5 


-5.5±0.5 


0.6±0.1 


1.7±0.1 


0.070 


B 


2135.9±0.1 


-4030.5±0.1 


3.3±0.2 


4.3±0.2 


0.065 


C 


1429.0±0.5 


-2889.9±0.5 


0.5±0.1 


0.5±0.1 


0.065 



2.3 Global VLBI 5 GHz observations 

In order to better resolve the north-east extension in image B and 
to determine the spectral index of component C, we undertook a 
global VLBI observation of B2108+213 at 5 GHz. The lens sys- 
tem was observed on 2006 February 17 with the Effelsberg, Jo- 
drell Bank Mk2, Westerbork, Medicina, Torun, Noto, Green Bank 
and the 10 VLBA antennas. The 14.5 h long observation was 
taken in both the right and left-hand circular polarizations, through 
4 IFs each with 8 MHz bandwidth and a 256 Mbs -1 bit rate. 
Here, the 5 min time cycle between the lens and the calibrator 



(B2103+213) was divided into 3 and 2 min scans, respectively 
for phase-referencing. The correlation was carried out at IIVBj 
where the data were divided into 16x0.5 MHz channels and time- 
averaged over 1-s intervals. The maximum field-of-view defined 
by bandwidth and time-averaged smearing is ~ 8 arcsec with this 
setup. Therefore, we obtained a single correlation of the data with 
the phase centre at the mid-point between lensed images A and B. 
The data were reduced using AIPS in a similar manner to the 



4 Joint Institute for VLBI in Europe 



© 2002 RAS, MNRAS OOO.rrHTTI 



Probing a massive radio galaxy with gravitational lensing 5 



Table 3. Positions and flux densities of the fitted Gaussian components for 
the 5 GHz data. Component A3 was not fitted with a two-dimensional Gaus- 
sian (see Section |2~3l for details). 



Comp. 


RA 


Dec 




Stotal 




(mas) 


(mas) 


(mJy beam - 1 ) 


(mJy) 


All 


0±0.1 


Q±0.1 


3.6±0.4 


3.8±0.5 


A12 


-u±o.i 


O.OrbO.l 


1.9±0.3 


2.6±0.3 


A21 


-3.5±0.1 


-0.6±0.1 


0.7±0.2 


1.2±0.2 


A22 


~4.7±0.1 


-0.5±0.1 


0.3±0.1 


0.3±0.1 


A23 


-6.0±0.1 


-1.3±0.1 


0.5±0.2 


1.4±0.2 


A3 


-9.8±0.5 


-3.9±0.5 


0.4±0.2 


0.4±0.2 


Bl 


2135.2±0.1 


-4030.0±0.1 


2.3±0.3 


2.9±0.4 


B2 


2137.0±0.3 


-4028.0±0.3 


0.5±0.2 


1.1±0.2 


B3 


2140.5±0.4 


-4025.2±0.4 


0.3±0.1 


0.3±0.1 


C 


1430.7±0.1 


-2888.1±0.1 


0.9±0.1 


0.9±0.1 



1.7 GHz VLB A observations. The Green Bank telescope had severe 
problems throughout the observing run, hence it was removed com- 
pletely from the dataset. All baselines with Torun were discarded 
because of a poor amplitude calibration. Throughout the data re- 
duction process all of the antennas were given equal weights to stop 
the large antennas (like the 100 m Effelsberg telescope) dominat- 
ing the uv-dataset. This avoids higher side-lobes and a deterioration 
of the image quality. The data for B2 108+2 13 were self-calibrated 
with a 3-min long solution interval, and CLEANed using IMAGR. 
The naturally weighted global VLBI maps of B2108+213 are pre- 
sented in the lower panel of Fig. [2] 

Image A shows the same core-jet structure observed at 
1.7 GHz and previously at 5 GHz (McKean et al. 2005). These new 
Global VLBI data clearly resolve image A into three main com- 
ponents Al, A2 and A3. Component Al is made up of two com- 
pact sub-components All and A12. Component A2 appears to be 
an extended knotty jet feature, which is further divided into three 
sub-components, A21, A22 and A23. The third component (A3) 
of image A is a faint jet feature which kinks towards the south. 
Components Al and A2 were fitted with two and three elliptical 
Gaussian models respectively, using Powell's minimization method 
jPressetal.ll 19921) in the image plane. Component A3 could not be 
well fitted with multiple Gaussian models. Therefore, the flux den- 
sity was measured by summing all of the surface brightness emis- 
sion within the 3-er boundary. The position of A3 was obtained 
from the surface brightness peak. Note that due to the different an- 
gular resolutions and frequency dependent structure, the compo- 
nents (Al, A2 and A3) at 1.7 GHz do not correspond to those ob- 
served at 5 GHz - this is simply our naming convention. The map 
of image B shows the expected core-jet structure to the north-east. 
We identify Bl as the core, and B2 and B3 as the counterparts to 
the jet features detected from the 5 GHz imaging of A. The emis- 
sion from image B was fitted with three elliptical Gaussian compo- 
nents, using the same minimization method as for image A. Note 
that the small feature to the south of component B2 was excluded. 
Component C was fitted with a compact elliptical Gaussian model 
(deconvolved FWHM is 0.8 mas) at 5 GHz. There is no evidence 
of any collimated jet emission towards the extended lobes detected 
in the 1.4 GHz MERLIN observations. 

The fitted positions and flux densities of the Gaussian model 
components are presented in Table[3] The errors for the flux densi- 
ties of A and B were calculated from the principle used in lFomalontl 
( 1999), and the errors on the positions were taken from the AIPS 
task JMFIT. 



3 MASS MODEL 

We now investigate mass models for the B2108+213 lens system. 
The current best model has two singular i sothermal spheres fo r 
galaxies Gl and G2, and an external shear (Mc Kean et alj|2005l) . 
This model reproduces the observed positions and flux densities of 
the lensed images, but has no degrees of freedom due to the limited 
number of observational constraints. Using the new high resolution 
VLBI observations presented here, we first test whether the image 
configuration can be better explained by a single lensing galaxy 
(Gl), or by two galaxies (Gl and G2). We then test whether com- 
ponent C can be a core lensed image. Finally, we investigate the 
influence of the environment on the mass model. The lens mod- 
els were gen erated using the publicly available code LENSMODEL 
dKeetorJl200ll) . 



3.1 Case 1: C as the active nucleus of lensing galaxy Gl 

The position of radio component C was used as the position for 
the lens galaxy Gl in this case. We note that the position obtained 
for component C from the VLBI imaging is not consistent with the 
position measured for Gl from the /- and //-band imaging (F814W 
and F160W, respectively) with the Hubble Space Telescope. This is 
almost certainly due to Gl having a complicated surface brightness 
profile, particularly in the infrared (this will be discussed further 
by McKean et al. in prep.). However, our position for component 
C is consistent with the V-band (F555W) imaging. Therefore, the 
position of G2 and the flux ratio between Gl and G2, which are 
used in the models, are taken from the V-band data (see Table 4 of 
iMcKean et al1l2005t) . 

As expected from lensing, we believe that the three main com- 
ponents of image A are the counterparts of the components in im- 
age B. However, image A is further resolved and shows several sub- 
components in Al and A2. The compact sub-component All is 
probably the core, hence it was used as the position for Al. As com- 
ponent A2 is much more extended, a flux-density weighted position 
was determined from the sub-components (Act — — 3.8±0.2 mas, 
AS = — 0.7±0.2 mas relative to All). For component A3 and all 
of the components of lensed image B, the positions from Table [3] 
were used. The flux densities for only Al (i.e. A11+A12) and Bl 
were used for the modelling. The total number of constraints to the 
lens models provided from the lensed images is 13. 

First, we consider the lens as a single galaxy Gl, with a singu- 
lar isothermal sphere (SIS) density profile, with an external shear. 
This model has 9 parameters (1 for the Einstein radius, 2 for the 
shear and position angle, and 6 parameters for the three source po- 
sitions) and 4 degrees of freedom (dof). The total xl of the model 
is 7.6, where the subscript refers to the dof, and the results are given 
in Table|4] This simple model fits the data well, and requires a large 
Einstein radius of 2. 1 8 arcsec and an external shear of 0.05 at a po- 
sition angle of 107.3 deg. 

We now test what contribution G2 has if it is included as part 
of the lens potential. Hence, we introduce galaxy G2 as a compan- 
ion lensing galaxy at the same redshift as Gl. A SIS centred on the 
optical pos ition of G2 was add ed to the SIS+shear model for Gl. 
Following McK ean et ail f2005), the mass ratio between Gl and G2 
was fixed to &G2/&G1 = 0.29±0.05, from the optical luminosities. 
Note that including G2 in this way does not change the number of 
dof. The total \\ of the model is 6.9 which is only marginally better 
than the single SIS+shear model used for Gl. The fitted parameters 
are given in Table|4] Including G2 lowers the Einstein radius of Gl 
to 1.71 arcsec and reduces the shear to 0.04. 
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Note that if the ratio of the masses is not fixed and the Einstein 
radius of G2 is left as a free parameter, then the Einstein radius of 
Gl is lowered further to 0.76 arcsec and the shear is reduced to 
0.01. Also, the Einstein radius of G2 increases to 1.5 arcsec. Here, 
the simplest model that can be fitted needs to have the main lensing 
galaxy positioned collinearly between the two lensed images. Since 
G2 is closest to this, it is taken as the main perturber and the shear 
is reduced. The reduced \ 2 of this model is 0.9, which is lower than 
for the case when the mass ratio is fixed (reduced \ 2 = 1-7) and for 
the single galaxy model (reduced \ 2 = 1.9). In Fig. [3] we show the 
X 2 plot of the Einstein radii of Gl and G2. There is a clear degen- 
eracy between the Einstein radii of Gl and G2. However, since Gl 
is clearly the dominant mass clump of the system as shown from 
the optical data, the prior from the optical luminosity is used for 
the SIS+SIS+shear model to break this degeneracy and to produce 
a physically plausible model. 

So far, only an external shear was used to account for any mass 
ellipticity in the model. It is clear that some ellipticity is required 
to account for the non-collinearity between the lensed images and 
lensing galaxy Gl. The ellipticity in the model could be due to Gl 
having an elliptical halo. Hence, the external shear was replaced 
with an elliptical mass distribution for Gl, while an SIS is fixed 
at the position of G2. The dof of this singular isothermal ellipsoid 
(SIE)+SIS model is 4, and the resulting reduced \ 2 is This 
model fits the image positions slightly better than previous mod- 
els without changing the fitted parameters significantly (see Table 
|4). The ellipticity of Gl is found to be 0.135 at a position angle of 
105.5 deg. The ellipticity of the surface brightness profile of Gl, as 
measured from the H ST data, is 0. 14 at a position angle of 57 deg 
dMcKean et alj2005h . The offset between the position angle of the 
halo and the light distribution of Gl suggests that the environment 
may also be affecting the shape of the lensing potential. However, 
since there are only two lensed images of B2108+213 it is not pos- 
sible to constrain both an elliptical potential for Gl and an external 
shear. Moreover, since ellipticity and shear (and also the presence 
of another lensing galaxy, G2) will produce similar observable ef- 
fects, the actual origin of the asymmetry between the lensed images 
and the lensing galaxy Gl is difficult to determine. 



3.2 Case 2: C as the core lensed image 

Extended mass distributions are expected to produce an ex- 
tra lensed image near the posit ion of the lensing galaxy 
jDver & Roede^U980l : lBurkdll98lh . However, the mass density 
distribution close to the centre of the lensing galaxy will affect the 
magnification of the core lensed image. For example, an isother- 
mal density profile will completely demagnify the core lensed im- 
age. Since galaxie s are known to have g lobal density profiles close 
to isothermal (e.g. Koopmans et al. 2006), searches for core lensed 
imag es tend to concent r ate on asymmetric double image systems 
(e.g. iBovce etalj |20bl: IZhang etal] [20071. that is, lens systems 
where the A:B flux ratio is greater than 10:1. This is because in 
these cases the magnification of the core lensed image is high- 
est and these systems offer the best possible chances of detection. 
The only example of a known galax y-scale lens system with a core 
lensed image is PMN J1632-0033 jWinn et al.ll2004h . which is an 
asymmetric double with a flux ratio of 15:1. For B2108+213, the 
flux ratio of image A to B is about 2:1, so it seems unlikely that 
component C could be a core lensed image. However, given the 
possibility that the lens potential may be made up of several galax- 
ies, the overall density profile might be shallower than isothermal. 



Table 4. The fitted parameters with the respective 1-<t uncertainties for the 
single lensing galaxy case (SIS+shear), and the two lensing galaxies cases 
(SIS+SIS+shear and SIE+SIS) with the mass ratio of G2 and Gl fixed. The 
angular separations relative to image Al and the Einstein radii are in mil- 
liarcsecond. The position angle of the shear and ellipticity is measured in 
degrees east of north. The superscript '*' refers to the values resulting from 
the fitted models. 



Name Fitted values 





SIS+shear 


SIS+SIS+shear 


SIE+SIS 


Gl (a 8) 


1431il 


1431il 


1431±1 




— 2888i 1 


7RXR+- 1 


2888i 1 


G2 (a, 8) 




1273±5 
-1786±5 


1273±5 
-1786±5 


Gl Einstein 


2177±16.0 


1711±38 


1736±38 


radius 








G2 Einstein 




498±39 


506±39 


radius 








External Shear 


0.050±0.007 


0.038±0.005 




Position angle 


107.3ijj 


10S.6±i« 




Ellipticity 






0.135±0.017 


Position angle 






105.5^2.1 


Source 1 (a, 8)* 


989.4 


1059.0 


1057.6 




-2110.3 


-2058.8 


-2021.6 


Source 2 (a, 8)* 


987.7 


1057.6 


1056.2 




-2110.1 


-2058.3 


-2021.1 


Source 3 (a, &)* 


985.4 


1055.7 


1054.4 




-2111.3 


-2059.1 


-2021.7 


Sba (Sbi/SaiT 


0.64 


0.63 


0.63 



Therefore, we now explore different mass models to establish the 
nature of component C. 

Throughout this section th e position of galaxie s Gl and G2 
are taken from the K-band data dMcKean et alj|2005h and the mass 
ratio is fixed. The positions and flux densities of A, B and C are 
taken from the global VLBI 5 GHz data (Table [3j. For the isother- 
mal models presented in Section [3~T1 no core lensed image is pro- 
duced. Therefore, we investigate mass distributions for Gl that de- 
viate from isothermal close to the centre of the lensing galaxy, that 
is, by inserting a core radius at the inner part of the halo. 

We use a Non-singular Isothermal Sphere (NIS) halo for Gl 
which has an isothermal mass profile with a flat density core. The 
second lensing galaxy G2 is too far away for its core properties to 
significantly affect the magnification of the central image. There- 
fore, we choose a SIS profile for galaxy G2. An external shear is 
also included in the model. This new model has three additional 
constraints from the position and flux density of component C and 
10 free parameters to fit, thereby giving 6 dof. The model fits the 
position of component C as a core lensed image well, but cannot fit 
the flux density, resulting in a very high total \a of 88.6. The best- 
fitting core radius is 0.1 mas (= 0.5 pc) and gives a flux density of 
10" 6 mjy for the core lensed image, which is still several orders of 
magnitude fainter than the observed flux density of component C. 
This model fails because a much higher core radius is required to 
fit the flux density of C, but this will be at the expense of fitting the 
position of component C and the positions and flux densities of the 
lensed images A and B. Therefore, it seems very unlikely that the 
third radio component is a core lensed image because the observed 
flux-density is just too high. 
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Figure 3. A degeneracy between the Einstein radii of Gl and G2 for the 
SIS+SIS+shear model is shown in the upper panel. Without any constraint 
on the masses of the galaxies, the best-fitting model requires G2 to be four 
times more massive than Gl, which contradicts the optical luminosities. In- 
troducing a constraint on the mass ratio between Gl and G2 breaks this de- 
generacy as shown in the lower panel. The intersection point of the double- 
dotted straight lines marks the minimum in both plots. 



3.3 Limits on the density profile of Gl 

For dark matter dominated structures like galaxy clusters, the inner 
density profile of the mass distribution tends to be quite shallow 
(i.e. 7 ~ 1- 1.5). This is based on the results from n umerical simu- 
lations (e.g. lNavarro et alll9 96; Moor e~et al J 19981) and from grav- 
itational lensing (e.g. lSand et al .1 12005). In this section, we attempt 
to place constraints on the density profile of the lensing halo, with 
the aim of determining whether it is consistent with an isothermal 
mass profile, or a shallower dark-matter dominated profile. To do 
this, we assume that the lensed image is not coincident with com- 
ponent C. We then derive an upper limit to the flux-density of a core 
lensed image using the 3-cr flux-density limit determined from the 
rms map noise (Sc ^ 150 ^iJy). 

We introduce a group halo with a power-law density profile, 




Figure 4. The relative magnification of image A and a possible core lensed 
image C as a function of density power-law slope for a single lensing galaxy 
(Gl) in the absence of a supermassive black hole. The solid (red) line is for 
a core radius of 0.005 pc. The dashed (green), dotted (blue) and dashed- 
dotted (cyan) lines are for core radii of 50, 150 and 250 pc, respectively. An 
upper limit to the relative magnification between images A and C is shown. 



external shear. The VLBI 5 GHz data for lensed images A and B 
are used to constrain the model. The Einstein radius, external shear 
and position angle are then optimized for density profiles with 7 
between 1.4 and 2.7, and for each value of the core radius. In Fig. 
[4] we show the relative magnification between lensed image A and 
a possible core lensed image C as a function of power-law density 
profile slope. We have also marked the upper limit to the relative 
flux ratio of a possible core lensed image and image A. First, it 
is clear that steeper density profiles result in a more demagnified 
core lensed image, as expected. Also, increasing the core radius 
from 0.005 to 250 pc, increases the magnification of the core lensed 
image. Interestingly, the upper limit to the relative magnification 
requires that the density profile of any group halo be steeper than 
~1.5. 

The presence of a radio-loud AGN implies that there is also a 
supermassive black hole at the centre of the lensing galaxy Gl. 
The presence of a black hole is expected to increase the num- 
ber of core lensed images to two, or possibly demagnify a core 
lense d image completely, depending on the mass of the blac k 



dope 

hole dMao. Witt & KoopmansfeOOlklRusin. Keeton & Winij2005h . 
Using the stellar velocity dispersion of the lensing galaxy Gl 
(a v = 360 kms -1 ; McKean et al. in preparation) and the known 
correlation between black hole mass and the stellar velocity disper- 
sion of the host galaxy, we find that the black hole associated with 
the AGN within Gl should have a mass of ~10 9 Mq. Including 
this black hole as a point mass at the centre of the Gl halo results 
in any core lensed image being completely suppressed. 

Our mass models which have Gl as the only lensing galaxy 
(SIS), or also include G2 as a companion lensing galaxy (SIS) to 
Gl (SIS or SIE), fitted the positions of the lensed images well, 
but failed to recover the flux ratio; the observed flux ratio is 
Sb I Sa = 0.45 at 5 GHz whereas the modelled flux ratio is ~0.63 
(see Table |4](. This may suggest that the overall density profile of 



par centred on the position of Gl. In addition, we incorporate 
a core radius of 0.005, 50, 150 and 250 pc. The model includes an 
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Figure 5. The \ 2 as a function of density power-law slope for the main lens- 
ing galaxy (Gl). The solid (red) line is without the inclusion of galaxy G2. 
The dashed (green), dotted (blue), dashed-dotted (cyan) and double-dotted 
(black) lines are for models that include G2 as an SIS with an Einstein ra- 
dius of 0.3, 0.5, 0.7 and 0.9 arcsec. 
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Figure 6. The relative magnification between images A and B as a function 
of density power-law slope for a single lensing galaxy (Gl). The solid (red) 
line is without the inclusion of galaxy G2. The dashed (green), dotted (blue), 
dashed-dotted (cyan) and double-dotted (black) lines are for models that 
include G2 as an SIS with an Einstein radius of 0.3, 0.5, 0.7 and 0.9 arcsec. 
The observed flux ratio at 5 GHz is shown with a conservative 20 per cent 
uncertainty in the flux densities of images A and B. 



the lens system differs from isothermal, possibly due to the sur- 
rounding environment. Therefore, we now investigate whether a 
different density profile is a better fit to the data. For this model, 
we again use a variable power-law density profile for Gl. As we 
are no longer using an isothermal mass profile for Gl, it is not pos- 
sible to fix the ratio of the Einstein radii of Gl and G2 from their 
optical luminosities via the Faber-Jackson relation. Therefore, we 
include G2 as a singular isothermal sphere with a fixed Einstein ra- 
dius of 0, 0.3, 0.5, 0.7 and 0.9 arcsec to test different scenarios. An 
external shear is also included. This model has 9 parameters and 
13 constraints. In Fig. [5] we show the model x 2 as a function of 
power-law density slope and for different Einstein radii for G2. It 
is clear that those models with a steeper density profile are a better 
fit to the data. However, increasing the Einstein radius of G2 also 
improves the fit, but not to the extent as to justify shallower mass 
profiles for Gl. Note that for the cases where the Einstein radius 
of G2 is 0.7 and 0.9 arcsec, shallow mass profiles for Gl (i.e. 7 ^ 
1.7 and 1.5, respectively) produce four lensed images due to the 
introduction of a new minimum and saddle point in the time-delay 
surface. Also, the flux density of image A is lower than that of im- 
age B. Therefore, these models can be ruled out. In Fig. [6] we show 
the predicted flux ratio of images A and B as a function of the den- 
sity profile of Gl and for the same set of Einstein radii for G2. The 
observed flux ratio at 5 GHz, with a 20 per cent error on the flux 
densities of images A and B, is also shown. It is clear that shallow 
density profiles are inconsistent with the flux ratio of A and B, if 
the halo is fixed at the position of Gl and assuming the flux densi- 
ties are not significantly affected due to variability and a time delay. 
Furthermore, models with shallower profiles (i.e. 7 ^ 1.7) require 
the flux density of B to be close to that of A, which is not consis- 
tent with the observed flux ratio. The best-fitting models require the 
density slope for Gl to be steeper than isothermal (7 = 2.45^o ig> 
for a single spherical mass model), and are fairly insensitive to the 
Einstein radius of G2. 



4 DISCUSSION 

The aims of this paper were to determine the nature of the third ra- 
dio component and to investigate the B2108+213 lensing potential, 
with particular emphasis on the contribution of the second lensing 
galaxy and the group environment. We now discuss to what extent 
these aims have been met. 



4.1 The lensing potential 

The new 1.7 and 5 GHz data presented here have found compact 
and extended emission from the two lensed images A and B which 
is consistent with gravitational lensing. We find that the surface 
brightnesses of the two lensed images are as expected for gravita- 
tional lensing. Furthermore, the series of three non-collinear com- 
ponents detected in the 5 GHz imaging of images A and B appear 
to have mirror symmetry. Hence, the expected parity reversal of the 
lensed images is also observed. In Fig. [7] we show the radio spec- 
tral energy distributions of images A and B. These spectra have 
been constru cted using the data p resented here and from the data 
presented bv lMcKean et al.l d2005l) . It is clear that both lensed im- 
ages have very similar flat radio spectra. We note that the MERLIN 
1.4 GHz flux-densities are much higher than the VLBA 1.7 GHz 
flux-densities for images A and B. This is either due to calibration 
errors in the MERLIN or VLBA data, or more likely, that extended 
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Figure 7. The radio spectra of lensed images A and B, and component C. 
The 1.4, 1.7 and 5 G Hz data points are from this paper. The rest of the data 
points are taken from McKean et alj j2005h . 



jet emission from the two lensed images has been resolved out by 
the high resolution 1.7 GHz imaging. 

We find that isothermal mass models which have only Gl as 
the lensing galaxy or which include the companion lensing galaxy 
G2, are a good fit to the data. However, the two lens galaxy model 
appears to be a slightly better fit. Both of these lens models fit the 
positions of the lensed images well, but do not fit the flux ratio 
well. This could be due to variability or a time delay in the radio- 
loud lensed images, or may indicate that a more complex mass 
model including additional group galaxies is required. Although 
B2 108+2 13 has not been monitored for variability there are 3 in- 
dependent observations carried out with MERLIN and the VLBA 
at 5 GHz to a low surface-brightness limit. As can be seen from 
Fig. [7J the flux densities from these observations are consistent to 
within their uncertainties. It seems likely that the lens environ- 
ment is playing a role in the image splitting of B2108+213. We 
find that dark matter dominated haloes, that is, those with shallow 
density profiles (7 ^ 1.5) seem not to be consistent with the ob- 
served flux ratios of the lensed images, or the non-detection of a 
core lensed image (in the absence of a supermassive black hole). 
In fact, density profiles that are steeper than isothermal are pre- 
ferred. Even though the density profil es of isolated galaxies appear 
to be consistent with isothermal (e.g. Koopmans et al. 200fJ), it is 
predicted that galaxies undergoing an interaction can have den- 
sity profiles that are steeper than isothermal for < 0.5 Gyr after 
the in itial interaction, before returnin g to the original isothermal 
state iDob ke. King & Fellhauei| [2007). Given the close proximity 
of galaxy G2, this scenario is certainly consistent with the lens sys- 
tem B2108+213. However, further observations will need to be car- 
ried out to confirm that the density profile of Gl is steeper than 
isothermal. In particular, measuring the redshift and stellar veloc- 
ity dispersion of the companion galaxy G2, coupled with the source 
redshift, will allow the inner density profile of the lensing mass dis- 



tribution to be determined. Also, if the background source yields a 
time-delay in the future, then the density profile can be found by as- 
suming the Hubble constant (e.g. |D obke & Kinj2 006; Auge r et all 



ing t 

l2007bl : lRead, Saha & Macci6ll2007t) . Since all of the mass models 
tested here require an external shear of ~4-5 per cent, it is likely 
that additional group/cluster galaxies are contributing to the overall 
potential. The impact of the group/cluster galaxies associated with 
Gl will be discussed in a follow-up paper. 

4.2 A radio-loud lensing galaxy 

It seems almost certainly the case, for several reasons, that the 
compact emission from radio component C is due to a radio-loud 
AGN hosted within the main lensing galaxy Gl. First, as can be 
clearly seen from Fig. [7] the radio spectrum of component C is ris- 
ing from 1.4 to 22.46 GHz, whereas the radio spectra of the lensed 
images A and B are flat. The spectral index of component C is 
Q 22.46 _ o.51±0.01. The spectral difference between component 
C and the two lensed images could be due to free-free absorption 
from the interstellar me dium of the lensin g galaxy (e.g. as in the 
case of PMN 1632-0033: Iwinn et al.ll2004l) . However, the effect of 
free-free absorption is known to be a strong function of frequency, 
and the spectrum of component C shows no curvature and is con- 
stantly rising from 1.4 to 22.46 GHz. Second, the radio emission 
is coincident with the optical position of the lensing galaxy Gl, as 
expected for an AGN. Third, although it is possible to fit a core 
lensed image at the position of component C, the flux density of 
component C is much too large to be a third lensed image. Note 
that this would be exacerbated if component C has been affected 
by free-free absorption. 

The classification of component C as the core of an AGN 
within Gl is also consistent with the extended lobe emission 
on either side of the lensing galaxy detected by MERLIN at 
1.4 GHz. Moreover, the diffuse lobed morphology of the ex- 
tended emissi on is as expected for a Fanaroff-Riley type I radio 
source (FRI; iFanaroff &RiIeyll 19741) . There are no highly colli- 
mated jets or hotspots which are typically seen in the more pow- 
erful FRII sources. The total 1.4 GHz rest-frame luminosity of 
the radio emission associated with Gl is L1.4 ~ 10 25 W Hz" 1 
(this assumes a spectral index of a = —0.8 for the lobe 
emission). Extrapolating this luminosity to 178 MHz gives 
Lo.i78 ~ 10 24 W Hz^ 1 sr" 1 , which is below t he FRI-FRII lu- 
mino sity divide of Lo.i78 ~ 10 25 WHz" 1 sr 1 (Fanaroff & Rilevl 
1 1974b . Finally, the absolute i?-band magnitude of the lensing galaxy 
Gl is Mr = —24. From the correlation between the rest-frame 
1.4 GHz luminosity and the absolute i?-ban d magnitude of the host 
galaxy (see fig. 1 of lLedlow & Owenl 19961) . which divides sources 
into FRI and FRII, we find that the Gl radio source is in the FRI 
region. Although component C appears to be unresolved in the 1 .4 
GHz MERLIN and 1.7 GHz VLBA imaging, the presence of the 
extended jet emission may explain why component C was s lightly 
extend ed in the MERLIN 5 GHz map presented bv lMcKean et al.l 
<2005h . 

There are now two gravitational lens systems from the CLASS 
survey with a known radio-loud lensing galaxy; the othe r is CLASS 
B2045+265 jFassnacht et alJl999l : lMcK ean et alj2007h . This gives 
a fraction of 1 1 ±8 per cent for lens galaxies with a radio-loud AGN 
from CLASS. This agrees closely with the percentage found from 
deep radio imaging of optically selected gravitationa l lens candi- 
dates from the Sloan Digital Sky Survey (~ 10 per centi lBovce et al.l 
120061) . 

It is possible to probe different lines-of-sight by searching 
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for differences in the properties of lensed images, which should 
be identical in the absence of variability. Since image separa- 
tions are typically only ~0.5-1.5 arcsec, it is the lensing galaxy 
which is mostly being probed. This technique has been most suc- 
cessful at optical wavelengths where the dust extinction along the 
lines-of-sight to each of the lensed images has been used to test 
galactic extinction laws in high redshift lens gala xies dFalco et all 
ll999l : IWucknitz et al.ll2003l ; lEliasdottir et al. 2006). At radio wave- 
lengths, free-free absorption of one (or more) of the lensed im- 
ages has probed the ISM of the lensing galaxy ( e.g. Iwinn et all 
120031 iMittal et alj|200l Iwiklind & Combed! 19961) . In the case of 
B2108+213, we have the unique situation of a radio jet from an 
AGN within the lensing galaxy passing in front of one of the lensed 
images (see Fig.Q}. Therefore, B2108+213 can be used to probe the 
composition of an FRI radio jet. As can be seen from Fig. [7] the 
flux ratio between images A and B is almost constant from 1.4 to 
8.46 GHU As such, there is no evidence of free-free absorption of 
the flux density of image B. Another possible propagation effect is 
depolarization/Faraday rotation of the emission from lensed image 
B as it passes through the magnetized plasma of the radio-jet. How- 
ever, the MERLIN data presented here were not deep enough to es- 
tablish whether or not image B has been depolarized. We do know 
that the background source is polarized by ~4 per cent at 1 .4 GHz 
from the data obtained for lensed image A. Therefore, further ra- 
dio imaging should be carried out to determine if the properties of 
image B have been affected by the radio jet. 



5 CONCLUSIONS 

We have presented new high resolution MERLIN and VLBI imag- 
ing of the gravitational lens system B2 108+2 13, whose wide image 
separation is consistent with a massive lensing galaxy or possibly 
a group of galaxies. The VLBI imaging at 1.7 and 5 GHz found 
extended emission in the lensed images whose surface brightness 
and parities are consistent with gravitational lensing. Using the new 
constraints provided from the two lensed images we have tested 
mass models for the B2108+213 lens potential. We find that the 
properties of the lensed images are consistent with either a sin- 
gle massive lensing galaxy, or a two galaxy lens model which ac- 
counts for a nearby companion to the main lensing galaxy. In these 
cases the companion galaxy G2 was always represented by a sin- 
gular isothermal sphere whereas the main galaxy was modelled 
as both an isothermal sphere and ellipsoid. We find that steeper 
than isothermal density profiles for the main lensing galaxy are 
a better fit to the data. Further models which include additional 
group galaxies will be tested when a spectroscopic survey of the 
B2 108+2 13 local environment has been completed. 

Emission from the third radio component, which is coincident 
with the optical position of the main lensing galaxy, was found 
for the first time at mas-scales at both 1.7 and 5 GHz. Further- 
more, MERLIN imaging at 1.4 GHz detected extended low surface 
brightness emission on either side of the third radio component. 
This jet emission has a morphology and luminosity which is con- 
sistent with an FR I type radio source. Attempts to model the radio 
core of the third radio component as a core lensed image failed be- 
cause it was not possible to fit the observed flux-density (to within 



5 The 22.46 GHz values have b een determined by as suming the flux ratio 
is unchanged at 22.46 GHz - see iMcKean et allfeOOSl) for details. 



several orders of magnitude) by using a core radius and/or a vari- 
able power-law density profile for the main lensing galaxy. There- 
fore, we are confident that the third radio component is due to emis- 
sion from an AGN embedded within the main lensing galaxy. 
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